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Abstract: - Partition functions Z{x) of statistical mechanics are generally approximated by integrals. The ap- 
proximation fails in small cavities or at very low temperature, when the ratio x between the energy quantum 
and thermal energy is larger or equal to unity. In addition, the exact calculation, which is based on number 
theoretical concepts, shows excess low frequency noise in thermodynamical quantities, that the continuous 
approximation fails to predict. It is first shown that Riemann zeta function is essentially the Mellin transform 
of the partition function Z{x) of the non degenerate (one dimensional) perfect gas. Inverting the transform 
leads to the conventional perfect gas law. The degeneracy has two aspects. One is related to the wave nature 
of particles: this is accounted for from quantum statistics, when the de Broglie wavelength exceeds the mean 
distance between particles. We emphasize here the second aspect which is related to the degeneracy of en- 
ergy levels. It is given by the number of solutions r^{p) of the three squares diophantine equation, a highly 
discontinuous arithmetical function. In the conventional approach the density of states is proportional to the 
square root of energy, that is r^^p) ~ 2iTp^^'^. We found that the exact density of states relates to the class 
number in the quadratic field Q(a/— p)- One finds 1// noise around the mean value. 

Key-Words: - Electronic circuits, number theory, quantum statistical physics, 1/f noise 



1 Introduction to Noise 
in Electrical Circuits 

Noise in electrical circuits is found in many 
forms, some of them have been well explained: ther- 
mal noise, shot noise, partition noise, burst noise... 
and one which is still subject to much debate due 
to its universality and the lack of a general accepted 
model: 1/f noise. 

Let us first review briefly our understanding of 
electrical thermal noise. Due to thermal agita- 
tion free electrons in a metallic conductor are mov- 
ing around continuously causing collisions with the 
atoms and a continuous exchange of energy between 
the modes. This was first investigated experimen- 
tally by Johnson ( I J and theoretically explained by 
Nyquist |2|. The noise in any circuit kept at uni- 
form temperature T can be described by a noise 
voltage(t'2)^/2 in series with a resistor R of the cir- 
cuit such that for a small frequency interval df 

^ = 4kRTp{f)df, (1) 

where p{f) = M(e'^//'=r _ is the Planck 
factor, k = 1.38 10"^^ J/K is Boltzmann's con- 
stant and h = 6.62 10"'^^ J.s is Planck's constant. 
For room temperature and not too high frequency 



hf /kT <^ 1, so that p{f) is equal to unity. These 
equations are verified from experiments. For exam- 
ple with a resistor of i? = 1 KQ, at room temper- 
ature T = 300 K, if one uses a filter of bandwith 
df = I MHz, the expected voltage fluctuation is of 
the order 4 fiV easily detectable on an oscilloscope. 

Our second example is shot noise in a diode. 
Such a noise is produced by random emission of 
electrons from the cathode to the anode. Schottky 
theorem [ 3J states that the mean current / undergoes 
a fluctuation 

^ = 2eldf, (2) 

where e = 1.6 10"^^ C is the charge of one elec- 
tron. Using a typical value 1 = 1 mA and the same 
bandwith as above the expected current fluctuation 
is on the order 18 nA. 

Our third exemple is partition noise which occurs 
any time a current is distributed randomly between 
two electrodes. Its net effect is an extra multiplica- 
tive factor in the relation for the shot noise current 
|4|. A modern variant of shot noise and partition 
noise is quantum partition noise. Eq. Q was used 
as an extra proof that the electrical charge is indi- 
visible and the measurement of current fluctuations 
allowed the determination of the unit e. Recently it 



has been argued that low dimensional systems like 
quantum wires or quantum dots may produce a frac- 
tional charge e| (p and q integers) associated to a 
quasi particle tunneling state. The corresponding 
fractional shot noise current has been observed [5 1. 

In most small size transistors a noise of the 
low frequency type is observed: burst noise also 
called random telegraph signal (or RTS noise). 
The noise typically consists of random discrete- 
switching events in the current flowing through a 
transistor. It was first seen in reverse-biased p-n 
junctions and bipolar transistors |6|. It was then re- 
lated to random single electron capture events into 
localized defects of very small transistors such as 
MOSFET's |7|. Assuming that to the probability 
of an electron making a transition from one state of 
amplitude to another state of amplitude I there is 
associated a switching time ti. The power spectral 
density is found to be of the Lorentzian type 

that is a flat spectrum up to the turn over frequency 
and a 1//^ dependance above. 

We now turn to 1// or flicker noise. Exam- 
ples above refer to a constant power spectral density 
Svif) = AkRT for thermal noise and Si{f) = 2el 
for shot noise. Flicker effect is the large ammount 
of noise generated in any solid state device (vacuum 
tube, diode, transistor...) at low frequencies which 
was discovered by Johnson 1 1 1. It is generally de- 
scribed phenomenologically as 

Svif) = (4) 

with a close to 2 and b close to 1. By integrating 
between the low frequency cut-off 1/r and the high 
frequency cut-off fc, we get 

^= f^^Svif)df^KVMfcr), (5) 
Ji/t 

that is the variability in flicker noise depends log- 
arithmically on the integration time r, as was ob- 
served experimentally |8|. Until now 1// noise in 
solid state systems eluded attempts to find its uni- 
versel origin. Fundamental experiments found 
flicker noise associated to thermal noise both in 
metals and semiconductors without a clear explana- 
tion on its microscopic origin. 

From many experiments it was found that the 
noise intensity is inversely proportional to the to- 
tal number of carriers A'^o in the sample that is 



K ~ 'y/No with 7 the Hooge factor fT(T|- Values 
of 7 in the range 10^'^ to 10~^ have been found. 
They were generally attributed to different scatter- 
ing mechanisms, by the lattice or the impurities, 
leading to mobility fluctuations of the electrons. 
This claims for a non linear origin of the l/f noise. 
Fine structures revealing the interaction of electrons 
with bulk and surface phonons were observed in dif- 
ferent solid-state physical systems by Mihaila |9J. 
On the theoretical side a quantum electrodynamical 
theory was developed by P. Handel |12| based on 
infrared divergent coupling of the electrons to the 
electromagnetic field in the scattering process. The 
basic result for the 7 parameter is as follows 



where 5w is the change of the velocity of the ac- 
celerated charge in the collision processes govern- 
ing the carrier mobility and the carrier diffusion 

constant, c = 3.10^ m/s is the velocity of light, 

2 

a = 2 hceo ~ 1/137 is the fine structure constant 
and eo = 8.85 10~^^ F/m is the permittivity of free 
space. Thirty years after its introduction the theory 
is still under debate llT3l.fl4l. 

In large devices such as MOSFET's l/f noise 
is often explained as arising from a summation of 
RTS's due to the defects in the insulator whenever a 
large range of time constants ti is involved. 

l/f noise is also observed in oscillators. One 
way to characterize it is from the Allan variance 

cr'yir) = ^S^, (7) 

where (5y(r) = 2/j+i(T) — yj(r) is the deviation be- 
tween two consecutive samples of index i and i + 1 
counted over an averaging time r and in = fi/fo 
refers to the ratio between the instantaneous fre- 
quency and the mean frequency of the oscillator. 
It can be shown [15| that l/f noise of the power 
spectral density Sy{f) = C / f, with C a constant 
related to the physical set-up, is associated to a con- 
stant value of the Allan deviation cTy(r) ~ 2 ln2 C. 
This so-called flicker floor is the limit to the stability 
of oscillators. Values in the range 10^^^ for quartz 
oscillators to 10^^^ for Hydrogen maser based os- 
cillators have been obtained fT6ll . 

We recently discovered a possible relationship 
between l/f frequency noise in oscillator measure- 
ments and prime number theory fFTI . lUSl . In this 
paper we pursue this quest by relating exact statis- 
tical mechanics of electrons of mass m in a box of 



size L to the quadratic field Q{^J—p) of negative 
discriminant, with p the number of energy quanta 
in units of 5E = h'^/SmL'^. We find numeri- 
cally a 1// noise of the fractional density of states 
about the average classical value. Sec. 2 shows 
that analytical number theory (and the Riemann zeta 
function) underlies the statistical mechanics of the 
perfect non degenerate gas. Sec. 3 shows that 
the degeneracy of energy levels needs the extended 
frame of algebraic number theory (and the field of 
quadratic forms). To let the text accessible to a 
large audience of readers we do not enter into the 
mathematical details that can be found elsewhere 

001,1111. 

2 The Non Degenerate Perfect Gas 

2.1 The classical approach 
to the perfect gas 

In statistical thermodynamics the partition function 
allows to establish all thermodynamical properties 
U91 . It is called generating function in the context 
of number theory. For a perfect classical gas the par- 
tition function per degree of freedom and per parti- 
cle is 



n=l 



(8) 



The energy levels are obtained by solving 
Schrodinger equation with free boundary con- 
ditions 

(fip 2mE , ^ 
= at X = and X = L. (9) 



The solution is 



^ = (_)i/2sin(-^), 



(10) 



with En = 6E n^, (n = 1, 2, • • •) so that the parti- 
tion function becomes 



Zq= ^ exp (-vrn^x), 



(11) 



n=l 



with X = 5E = j .Using the first order 
Mac-Laurin expansion and if x is small 

Zo{x) ^ ^erf(n_V^) - ^ - (12) 

For the three dimensional gas of N particles we 
need the indistinguishability parameter A^! so that 



h2 



z = \. 

m 

The gas pressure is defined from 



(13) 



P{V,T) = 3NkT{ 



dlnZ^ 
dV 



)t.n withV = L'^, 



(14) 



The law of perfect gas follows 

PV = RT withR = kNA, (15) 

where Na is the Avogadro constant. 

Others useful state equations are derived such as 
the mean energy per mole of the gas 



and the free energy F and the entropy S as 

F = -kT\nZ. 



(16) 



(17) 



The Sackur-Tetrode equation follows 



S = S{E, V, N) 



Nk 



| + ln(^) + |ln(^f; 



(18) 



Thanks to the inclusion of A^! term in the definition 
of partition function the entropy equation satisfies 
the property of extensivity. 

2.2 The classical law of the perfect gas 
and the Riemann zeta function 

The partition function Zq appears in Riemann's pio- 
neering work on prime number theory. It starts from 
the integral representation of gamma function 

r(s) = / exp(-x)x'~^d2; with K(s) > 0, 
Jo 

(19) 

which generalizes the factorial function into the 
complex plane s. In particular r(s + 1) = sT{s) 
and as result T{n + 1) = n! if n is an integer. 
Riemann's zeta function is defined as 

C(s) = E— with 3f?(s) > 1. (20) 



n=l 



n 



Substituting vrn^x to x in the integral (IT9l and sum- 
ming for all integers n one gets 

/•CO 

7r"'r(s)C(2s) = / Zoix)x'~^dx, (21) 
Jo 



which is the desired relation between (^{s) and 
the Melhn transform of Zq{x). Here we assumed 
iT'ma.x — > oo in the definition of the partition func- 
tion. The expression 



as) 



-/2r(./2)C(5), 



(22) 



is the so-called completed Riemann zeta function. 
This is because ^(s) admits an analytic continua- 
tion to the whole complex plane s except at the 
simple poles s = and s = 1 of residus —1 and 
1 respectively. It satisfies the functional relation 

The approximation of Zq{x) may be obtained by 
taking the inverse Mellin transform and by applying 
the residue theorem to the two poles of C(2s) at s = 
and s = 1/2 so that 



Zoix) ~ 



1 

2' 



(23) 



which is similar to the previous Mac-Laurin approx- 
imation iTH . The step from the discrete to the con- 
tinuous is condensed into the contribution of poles 
of the Riemann zeta function. 

The remaining poles at s = —m (m integer > 1) 
of r(s) are absorbed by the trivial zeros at s = —I 
(I a positive integer) of C{2s). The unproven Rie- 
mann hypothesis 1 17 1 is that all remaining zeros are 
located on the critical line s = 1/2. The modulat- 
ing action of the zeros is expected to be significant 
within a box at the nanoscopic scale x = ^ > 1. 



3 The Degenerate Perfect Gas 

3.1 First form of quantum degeneracy 
due to wave-particle aspects 

In the quantum gas, the first form of degeneracy is 
due to the indiscemability of (wave)particles of size 
Ath (the de Broglie wavelength) larger than the typ- 
ical spacing d. The degeneracy parameter is 1 19] 



Ath h/{mv) 27rkT 
^ = -^' withv=(-— ) 



1/2 



(24) 



(For conduction electrons in silver at room temper- 
ature the ratio is about 40). 

The partition function for electrons is taken to 
satisfy Fermi-Dirac statistics with chemical poten- 
tial a 

InZ = ^ In [1 + exp(-a - +iVa. (25) 



In the conventional approach of statistical mechan- 
ics, the discretness of quantum states is neglected 
and one goes from the discrete to the continuous by 
taking the number of states ^q{E) as |, times the 
volume of the Fermi sphere of radius Rp, times 2 
(for the two spin states of the electron). The density 
of states ujq{E) per unit volume V = follows 
from 

ME) - y^r^ = M-f^) ' El . (26) 

3.2 Second quantum aspect 

due to the degeneracy of energy levels 

Here the number of possible quantum states is cal- 
culated exactly from number theory. 



n{E) 



E 



hni<-R| 



1 with R| = (27) 



where ni, n2 and 723 are positive integers, Ep is 
the Fermi energy and 6E is the energy quantum as 
above. In quantum dots the energy is quantized 

in units of 6E: this may be observed when the 

2 

Coulomb blockade energy Eb = fj is such that 
Eb — 5E S> kT, that is at very low temperatures 
or for small samples with large capacity C. The par- 
tition function is 



Pmax 



E 

^5pexp(--^) withEp 

p=i 



p, (28) 



where the degeneracy parameter Qp is the number of 
states of energy Ep. 
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Figure 1: Relative fluctuations e(p) of the density of 
states versus energy p (We took 5p = 100). 

The number of quantum states gp also equals |, 
times 2, times the number rs (p) of solutions of the 
three squares diophantine equation 



2 , 2 , 2 
ni + + 



p with p < Rp, 



(29) 



and rii, n2, ,p integers. This is a central prob- 
10 b — — — I 
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Figure 2: Power spectral density of the fluctuations 

lem of number theory. Gauss proved that for p > 4 
and squarefree 

r 24/i(-p) forp = 3 (mods), 
i^sip) = { 12/i(— 4p) for p = 1, 2, 5 or 6 (mod 8), 
[ for p = 7 (mod 8), 

(30) 

where h{—p) is the class number in the field 
Q{s/—p) (that is the number of primitive quadratic 
forms with given negative discriminant -p). The 
general case is r3(p) = Y.d?\p Rsi^)' where Rs{p) 
holds for square free values of p. More explicit so- 
lutions can be found ED- They are available in the 
Mathematica package fTT\. 

The class number h{—p), and thus r^{p) is a 
highly discontinuous function. We defined the (ex- 
act) density of states as 

'^(p) = iEn=r'"^3(n) 

Relative fluctuations of the density of states 
around the mean value LOo{p) were calculated as 

<p) = ^^^^w^ ^'^^ Pig- ^ 

suit the power spectral density of low frequency 
fluctuations versus Fourier frequency / was found 
as FFT{f) ~ A|2, with a close to 1 and A ~ 0.1, 
that is a 1// spectrum as found in Fig. 2. 

Physically one can take 5p = [|^], with [ ] the 
integer part and 5E the granularity parameter, since 
kT is the width of the Fermi surface at finite tem- 
perature T. One can conclude that the 1// spectrum 
arising in the arithmetical approach connects to the 
different regimes one gets for the electron gas in 
quantum dots, versus the typical size L of the dot, 
the temperature T and the effective mass m. 



Due to the low frequency noise in the density of 
states all thermodynamical quantities are perturbed, 
as well as the conductivity of electronic devices. 
This will be developed in a future paper. 

4 Concluding Remarks 

Thermodynamical quantities attached to the quan- 
tum description of an electronic gas in a cubic box 
are not smooth versus energy or time, contrary to the 
classical viewpoint: this is because the summatory 
function of an arithmetical function such as r3(n) is 
discontinuous. The fractional deviation with respect 
to the asymptotic average very often leads to a low 
frequency power spectrum of the fractal type 1//'', 
h close to 1. Number theory thus is a good candidate 
to explain many low frequency noises encountered 
in solid state physics and engineering as well as in 
others fields. 
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